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A regularity theorem for solutions of the spherically 
symmetric Vlasov-Einstein system 
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Abstract 



We show that if a solution of the spherically symmetric Vlasov- 
Einstein system develops a singularity at all then the first singularity 
has to appear at the center of symmetry. The main tool is an estimate 
which shows that a solution is global if all the matter remains away 
' from the center of symmetry. 

o ■ 

O ■ 1 Introduction 

m ■ 



This paper is concerned with the long-time behaviour of solutions of the 
spherically symmetric Vlasov-Einstein system. In [4] a continuation crite- 
rion was obtained for solutions of this system, and it was shown that for 



5-H . 

small initial data the corresponding solution exists globally in time. In the 
fohowing we investigate what happens for large initial data. In the coordi- 
. ■ nates used in [4] the equations to be solved are as follows: 

X 

c5 ■ dtf + e>^-^^I=-d,J-(—X + e>^-^Vl + ^fi')--d.J = 0, (1.1) 

e^2^(2rA'-l)-M = 87rrV, (1-2) 
e~^^(2r/i'+l)-l = 87rr2p, (1.3) 



pit,x) = / VT+^fit,x,v)dv, (1.4) 
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Here x and v belong to IR^, r:=|3;|, x-v denotes the usual inner product of 
vectors in ]R^, and v"^ :=v-v. The distribution function / is assumed to be 
invariant under simultaneous rotations of x and v, and hence p and p can 
be regarded as functions of t and r. Spherically symmetric functions of t 
and X will be identified with functions of t and r whenever it is convenient. 
In particular A and p are regarded as functions of t and r, and the dot and 
prime denote derivatives with respect to t and r respectively. It is assumed 
that f(t) has compact support for each fixed t. We are interested in regular 
asymptotically flat solutions which leads to the boundary conditions that 

X(t,0) = 0, lira ij(t,r) = 0, (1.6) 

r — 5-00 

for each fixed t. 

The main aim of this paper is to show that if singularities ever develop 
in solutions of the system (|l.l| )-(p^) with the above boundary conditions 
then the first singularity must be at the center of symmetry. In order to do 
this we consider solutions of ( |1.1D -( |L5| ) on a certain kind of exterior region 
with different boundary conditions and prove that for the modified problem 
there exists a global in time solution for any initial data. One of the most 
interesting implications of the results of is that the naked singularities 
in spherically symmetric solutions of the Einstein equations coupled to dust 
can be cured by passing to a slightly different matter model, namely that 
described by the Vlasov equation, in the case of small initial data. The 
results of this paper strengthen this conclusion to say that shell-crossing 
singularities are completely eliminated. 

The following remarks put our results in context. For the Vlasov- Poisson 
system, which is the non-relativistic analogue of the Vlasov-Einstein system, 
it is known that global existence holds for boundary conditions which are 
the analogue of the requirement of asymptotic flatness in the relativistic case 
[§, m, P and also in a cosmological setting . No symmetry assumptions are 
necessary. For the relativistic Vlasov-Poisson system with an attractive force 
spherically symmetric solutions with negative energy develop singularities in 
finite time Q. It is easy to show that in these solutions the first singularity 
occurs at the center of symmetry. On the other hand it was also shown 
in Q that spherically symmetric solutions of the relativistic Vlasov-Poisson 
system with a repulsive force never develop singularities. The latter are 
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in one-to-one correspondence with spherically symmetric solutions of the 
Vlasov-Maxwell system. 

The paper is organized as follows. Sect. 2 contains the main estimates 
together with a proof that they imply global existence in the case that all the 
matter remains away from the center. In Sect. 3 a local existence theorem 
and continuation criterion for the exterior problem are proved. It is then 
shown that the estimates of Sect. 2 imply a global existence theorem for the 
exterior problem. Finally, in Sect. 4, these results are combined to give the 
main theorem. 



2 The restricted regularity theorem 

The goal of this section is to show that a solution may be extended as long as 
/ vanishes in a neighborhood of the center of symmetry. Consider an initial 
datum / >0 which is spherically symmetric, C , compactly supported, and 
satisfies 

f j Vl + v'^f{x,v)dvdx<r/2 (2.1) 

J\x\<r J 

for all r >0. By Theorem 3.1 of Q a regular solution {f,X,fi) of the system 
( |l.l[ )-(p!?^) with boundary conditions ( |1.6D and initial datum / exists on 
[0,r[xIR'^ for some r>0. 

o 

Theorem 2.1 Let f and T he as above (T finite). Assume there exists e > 
such that 

f{t,x,v)=OifO<t<Tand\x\<e. (2.2) 
Then (/, A,/i) extends to a regular solution on [0,T'[ for some T' >T. 

Define 



P 



(t) ■.= snp{\v\ : ix,v) G supp/(t)}, (2.3) 



then Theorem 3.2 of states that Theorem |2.1| above follows once P{t) is 
shown to be bounded on [0,r[. We need a few other facts from Q: 

e~^^ = l-2m/r (2.4) 

where 

m{t,r):=A-K s^p{t,s)ds (2.5) 
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(equations (2.11) and (2.12) of [§). Also 

A = -47rre''+^j (2.6) 

where 

j(t,r):=J^f{t,x,v)dv (2.7) 

(equations (3.37) and (3.38) of [§). 
The fohowing notation wih be used: 

ii:=|v|, w:=r~^x - v, F ■.= \x Av\'^ = r'^u'^ — {x ■v)'^ . 

Differentiation along a characteristic of the Vlasov equation is denoted by 
Dt, so 

and 

\ r J r 

It follows that 

DtF = 0, (2.8) 
Ar = e^-^^=^, (2.9) 

and 

Dtw = —r^'^{Dtr)x-v + r~^{Dtx) ■v + r~^x- {Dtv). 
Substitution for Dtr,Dtx, and DfV and simplification yields 

Dtw = ^=ef'-^-wX-ef'-^Vl + ^f^'. (2.10) 

The letter C will denote a generic constant which changes from line to line, 

o 

and may depend only on f,e, and T. 

Now we make a few preliminary estimates. The values of / are conserved 
along characteristics so 

0</<sup/ = C7. 



Also we claim that 



J p{t,x)dx = J p{0,x)dx = C (2.11) 
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To show this multiply ( |1.1D by Vl + v^ and integrate in v, which yields (after 
simplification) 

= dtp + ef'-^div (^j fvdv')+{p+p)X + 2jef'-^n' 

= dtp + div {e^'^J fv dv) + {p+p)\+ je^-^ {p' + A') • 



Now substituting ( ^ ) and (1^), (O) this becomes 



Q = dtp + di-v(e^'-^ j fvdv) (2.12) 
and pup follows. Also by 

0<m{t,r)< J p{t,x)dx<C, r>0. (2.13) 



It follows from (O), (O) that 



and from (1.6) and ( |2.4| ) that 

lim (u + A) = 0, 

SO 

/i-A</i + A<0 

and 

gM-A<gM+A<i^ (2.14) 

Note that 

0<F<C 

on the support of /, so 

u^ = yj2^^<nj^ + ^ = ^2^C. (2.15) 

Hence we will focus on w. Define 

Pi{t) :=inf|ti; : 3x,v with f{t,x,v) and w = r^^x-v^ (2-16) 

and 

Ps{t) := sup|it; : 3a;, f with f{t,x,v) / and w = r~^x-v^, (2-17) 



then if -Pj(t) and Ps{t) are bounded, it follows that P{t) is bounded. Also 
note that 

measure \^v:{x,v)£suppf{t)^<TTCe''^(^Ps{t)-Pi{t)y \x\>e. (2.18) 

Next we focus on the characteristic equation for w. Note first that by 
(PD and dU) 



^'=-( e^\87Tr^p+ l)-l)= -e^^ ( 87rr2p+ 1 - 1 + — ) 
= e^'^(r~^m + 47rrp). 
Using this and ( |2^ ) in ( |2.10| ) yields 



F 



Dtw = e^~^-r~^\/l + n2e^+^m + 47rre^+^(wj - \/l + v^p) 

(2.19) 

By (2.2) and (2.14) we may bound the first term of ( 2.19| ) by 



< ^ e^-^ < e-^C = C 



on the support of /. By (|1|), (|2l^ , and (|2.13D we may bound the 

second term of ( p.l9| ) by 



< r-Vl + u2e^+^m < e^VC + w^c < Cv^C + w^ 

on the support of /. Hence, ( p. 191 ) becomes 

-CVc + w^+inref'+^iwj-VT+v^p) < Dtw (2.20) 

< C + 47rre^+^(u;j- VT+n2p). 

On the support of / 

< 47^re^'+^ < C 
so we must consider the quantity wj — ^/l + vflp. Let us denote 

— |~| 771 I . ~|2 2-2 / ~\2 

w:=r x-v, u:=\v\, ■.= \xAv\ =r u —[x-v) . 

Then 



/^2 
f{t,x,v)(ww — \/l + u'^—^ „ )dv. (2-21) 



The next step is to use and (Q) to derive an a upper bound for 

w (on the support of /), and hence for Ps{t). This may be done without a 
bound on Pi{t). Then ( 2.2C| ), ( p. 21 ), and the bound for Ps{t) will be used to 
derive a lower bound for w, and hence for Pi{t). Then by ( p.lSD and (|2.3D a 
bound for P{t) follows, and the solution may be extended. 

To bound Ps{t) suppose 

Ps{t)>Q 

and consider w (in supp/) with 

w>Q. 



For tt; < we have 



For > we have 



WW - 



WW ■ 



w 



w 



<0. 



w w'^ {1 + u'^) — w'^ {1 + u^) 



w w'^{l + Fr-'^)-w^{l + Fr 
wV^ + u'^ + wV^ + u'^ 



Note that in the last step a term of ''^w w canceled, which is crucial. Hence 



WW - 



w'^ 



< 



w 



w'^il + Ce' 



<C- 



WW 



1 + w 



7,2 ■ 



Now using the above and ( 2.1^ ) we have 



J f(ww- 



w 



dv < 



fCw 



w 



0<W<Ps{t) 

Psit) 



l + w 



r,2 



dv 



Cw 



w 



1 + w'^ 



dw 



= Cwln{l + P^{t))<CPs{t)ln{l + P^{t)), 

and by 1^1^ and (|]2^) 

Dtw<C + CPs{t)ln{l + P^{t)), (2.22) 

for > in the support of /. Denote the values of w along a characteristic 
by w{t) and let 



to:=inf{T>0:w(s)>0 for sG]r,t[}. 
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Then either to = or w(to) = and in either case 

w{to)<C. 

Hence by ( pl 

W{t) < w{to)+j'{C + CPs{T)Hl + P^{T)))dT 

<C + C f Ps{T)\n{l + P'^{T))dT. 
Jto 

Defining 

P,(r):=max{0,P,(r)} 

we may write 

w{t)<C + C I Ps{T)\n{l + P]{T))dT 

Jo 

and hence 

Ps(.t)<C + C f'Ps{T)\n{l + Pl{T))dT. 
Jo 

We assumed that Ps{t) > 0, but note that the last inequality is valid in all 
cases. It now follows that 

Ps{t)<Ps{t)<exp{e^^)<C (2.23) 

onte[0,T[. 

To bound Pi{t) from below suppose 

Pi{t)<0 

and consider w (in support /) with 

Pi{t)<w<0. 

For 'w<0 we have 

w w'^{l + Fr''^)-w'^{l + Fr-^) 



WW — u 



_ \w\ w'^{l + Fr''^)-w'^{l + Fr~'^) 
~ VlT V? \w\\/l + v? + \w\\/l + v? 
\w\ (-t()2)(l + Fr~2) 



> 
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^ -w^{l + Ce- 



'w\ 



> -C 



For w<0<w<Ps{t) we have (using (|2.23| )) 



w'^ „ , n. — ^ nP' 



WW - 



Vl + v?—== > Ps{t)w-\/l + w^ + Fr 



> Cw-\w\Vl + Ce-^ + w^ 

> Cw-Ps{t)VC + w^ 



w 



2 



> -cVc+w'^. 



Hence (using (|2.1g| ) as before) 

r,2 



f(ww — \/l + u'^ , \dv 



> / H-c p \dv^ f f(-cVc+wAdv 

Jw<0 ^ Vl + ^y^ Jw>0 ^ ' 

n rO rPs{t) 

> , 7rCe~^ / \w\dw-CVC + w^7rCe'^ 

VI + W^ J-Piit) Jo 



dw 



1 



CP^ (t) ^ - CP, (t) VC + w^. 



Vl + w 



Now by (II21I ), ( ^) , and (|2^ ) 

Aw^ > -C y/C + w2 - C (t) — - CPs it)VC + w^ 

Vl + w^ 

> -cVc+^-cpf{t)^l= 



^/l + w^' 

Since we have assumed > > Pj(t), it is convenient to write this as 
Dtiw"^) = 2wDtw 

i-w) 



< C{-w)VC + w^ + CP^{t) 

<cm)\^/c+pm+cpHt) 

< C + CPf{t). 



w" 



As before define 



then 



so 



ti:=mf|T>0:u;(s)<0 for se]r,t[} 
0>u;(ti)> -C 

w\t)<C+ l\c + CPf{T))dT 
Jti 

<C + C f Pf{T)dT. 

Jo 



It follows that 



Pf{t)<C + c( Pf{T)dT 

Jo 

if Pi{t) < 0. But if Pi{t) > then 

0<P^{t)<P,{t)<C 



(2.24) 



so (2.24) holds in this case, too. Now by Gronwall's inequality it follows 
that 

Pi{t)<e''<C 



on [0,T[. Finally a bound for P{t) follows from (|2.15D and (|2.3D , and the 
proof is complete. 



3 The exterior problem 

If ri and T are positive real numbers, define the exterior region 

W{T,ri):={{t,r):0<t<T,r>ri+t}. 

In this section the initial value problem for ( |l.l[ )-( pr^ ) will be studied on a 
region of this kind. Consider an initial datum f{x,v) defined on the region 
|a;|>ri which is non-negative, compactly supported, C^, and spherically 
symmetric. The first of the boundary conditions ( |1.6| ) cannot be used in the 
case of an exterior region, and so it will be replaced as follows. Let rrioo be 
any number greater than Att J^r'^p{0,r)dr. For any solution of (|1.1| )~( |L5|) 
on W{T,ri) define 

roc 

m(t,r) := moo — 47r / s'^p{t,s)ds. (3-1) 

Jr 
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Provided moo satisfies the above inequality, the quantity m(0,r) is every- 
where positive. The replacement for (|1.6|) is: 



"2^(*''') = l-2m(t,r)/r, lim /i(t,r) = 0. (3.2) 



The first of these conditions is a combination of ( |1.2| ) with a choice of bound- 
ary condition. Note that if a solution of the original problem with boundary 
conditions ( pT^ ) is restricted to W{T,ri) then it will satisfy ( |3.2D provided 
nioo is chosen to be equal to the ADM mass of the solution on the full space. 
Just as in the local existence theorem in [4], a further restriction must be 
imposed on the initial datum. In this case it reads 



|a;|>r 



\/l + v^f{x,v)dvdx <r/2, r>ri. (3.3) 



This is of course necessary if ( |3.2D is to hold on the initial hypersurface. 
The nature of the solutions to be constructed is encoded in the following 
definition. 

Definition A solution (/, A,/i) of ( |l.lD -(pr5|) on a region of the form W{T,ri) 
is called regular if 

(i) / is non-negative, spherically symmetric, and C^, and f{t) has com- 
pact support for each t£ [0,T[ 

(ii) A>0, and A, /u. A', and fi' are C^. 

A local existence theorem can now be stated. 

o 

Theorem 3.1 Let moo>0 be a fixed real number. Let f>0 be a spheri- 
cally symmetric function on the region \x\ > ri which is and has compact 
support. Suppose that i\3.!^ ) holds for all r>ri and that 



\J 1 + f{x,v)dvdx < moo. (3-4) 

/|a;|>ri J 

Then there exists a unique regular spherically symmetric solution of ( [i. 

o 

(1.5) on a region W{T,ri) with /(0) = / and satisfying (S.i. ). 

Proof : This is similar in outline to the proof of Theorem 3.1 of [4] and 
thus will only be sketched, with the differences compared to that proof 
being treated in more detail. Define Ao(t,r) and fj,Q{t,r) to be zero. If 
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A„ and Hn are defined on the region W{Tn,ri) then fn is defined to be 
the solution of the Vlasov equation with A and u replaced by and Un 

o 

respectively and initial datum /. In order that this solution be uniquely 
defined it is necessary to know that no characteristic can enter a region of the 
form W{T,ri) except through the initial hypersurface which is guaranteed if 
An > and < on the region of interest. This will be proved by induction. 
If /„ is given then A„+i and fin+i are defined to be the solutions of the 
field equations with p„ and Pn constructed from /„ rather than /. The 
quantities A„+i and fin+i are defined on the maximal region W{Tn+i,ri) 
where < m„(t,r) < r/2 so that A„+i can be defined by the first equation 



in (3^) and is positive on VF(T„+i,ri); note that < m„(t,r) < r/2 for r 

o 

large and < 771^(0, r) < r/2, r > ri by assumption on / so that T„+i > by 
continuity. It can be shown straightforwardly that the iterates (A„,/in)/n) 
are well-defined and regular for all n. The most significant difficulty in 
proving the corresponding statement in [4] was checking the differentiability 
of various quantities at the center of symmetry, and in the exterior problem 
the center of symmetry is excluded. The next step is to show that there 
exists some T > so that T„ > T for all n and that the quantities A„ , A„ , 
and /u^ are uniformly bounded in n on the region W{T,ri). Let Ln{t,r) := 
r-i(l-e-2-^"(*'''))-i. For te[0,Tn[ define 

Pn{t) := sup{|t;| : {x,v) £suppfn{t)}, 
Qn{t) := ||e2^"W||oo + ||L, 



Now it is possible to carry out the same sequence of estimates as in [4] 
to get a differential inequality for Pn{t) ■.= niaxQ<k<nPkit) and Qn{t)'-= 
maxo<fe<nQfe(i) which is independent of n, and this gives the desired result. 
The one point which is significantly different from what was done in [4] is the 
estimate for A^. There a partial integration in r must be carried out, and 
in the exterior problem the limits of integration are changed; nevertheless 
the basic idea goes through. The remainder of the proof is almost identical 
to the proof of Theorem 3.1 of [4]. It is possible to bound /i^ and A^ on 
the region W{T,ri) and then show that the sequence {Xn,IJ'n,fn) converges 
uniformly to a regular solution of ( |l.lD -(pr5|) on that region. Moreover this 
solution is unique. □ 

Next a continuation criterion will be derived. By the maximal interval 
of existence for the exterior problem is meant the largest region W{T,ri) on 
which a solution exists with given initial data and the parameters ri and 
nioo fixed. 
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Theorem 3.2 Let (/, A,^) be a regular solution of the reduced system 

. o 

on W{T,ri) with compactly supported initial datum f. If T< 00 and 
W{T,ri) is the maximal interval of existence then P is unbounded. 

Proof : Note first that, just as in the case of the problem in the whole space, 
the reduced equations ( |l.lD - (ll.5|) imply that the additional field equation 
(|2.6|) is satisfied. The conservation law (|2.12|) can be rewritten in the form 



^t{ry) + ^r{e^'-^r^j) = {). (3.5) 

Integrating (^]^) in space shows that the minimum value of m at time t 
(which occurs at the inner boundary of W{T,ri)) is not less than at t = 0. 
Hence L is bounded on the whole region. Suppose now that P is bounded. 



From (2.6) it follows that A is bounded, and if T is finite this gives an upper 
bound for A. Combining this with the lower bound for A already obtained 
shows that Q is bounded. This means that the quantities which influence the 
size of the interval of existence are all bounded on VK(r,ri), and it follows 
that the solution is extendible. □ 



Combining Theorem |3.2| with the estimates of Sect. 2 gives a proof that 
the solution of the exterior problem corresponding to an initial datum of 
the kind assumed in Theorem 13.11 exists globally in time i. e. that T can be 



chosen to be infinity in the conclusion of Theorem 3T. For we know that it 
suffices to bound P, and the estimates bound the momentum v along any 
characteristic on which r is bounded below. Moreover the inequality r > ri 
holds on Ty(oo,ri). 

4 The regularity theorem 

This section is concerned with the initial value problem for (|l.l|)- (|l.5| ) on 
the whole space with boundary conditions (|l7 



Theorem 4.1 Let (/, A,/i) be a regular solution of the reduced system 



(1.5) on a time interval [0,T[. Suppose that there exists an open neighbor- 



hood U of the point {T,0) such that 

sup{|?;| : (t,a;,f) Gsupp/n(C/ xIR^)} <oo. (4.1) 
Then (/, A,;u) extends to a regular solution on [0,T'[ for some T' >T. 
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Proof : Suppose that the condition ( ^ ) is satisfied. Since the equations 
are invariant under time translations it can be assumed without loss of 
generality that T is as small as desired. Choosing T sufficiently small ensures 
that U contains all points with {T - tf + < 4T'^ and 0<t<r. Now let 
n<T. Then [0,T[xM^ cULiW{T,n). Let / be the restriction of / to the 

o 

hypersurface i = 0. Restricting / to the region > ri gives an initial datum 

o 

for the exterior problem on W{T,ri). Let nice be the ADM mass of /. 
There are now two cases to be considered, according to whether / vanishes 
in a neighborhood of the point (T,0) or not. If it does then by doing a 
time translation if necessary it can be arranged that the matter stays away 
from the center on the whole interval [0,r[ so that the results of Sect. 2 are 
applicable. It can be concluded that the solution extends to a larger time 
interval in this case. If / does not vanish in a neighborhood of (T,0) then 
by doing a time translation it can be arranged that ( |3.4| ) holds. In this case 
the results of the previous section show that there exists a global solution 
on I^(oo,ri) satisfying (|3.2|). The extended solution must agree on W{T,ri) 
with the solution we started with. Thus a finite upper bound is obtained for 
\v\ on the part of the support of / over W(T,ri). But by assumption (4.1 



we already have a bound of this type on the remainder of the support of /. 
Hence 

sup{|f I : {t,x,v) G supp/} < oo. 

Applying Theorem 3.2 of [4] now shows that the solution is extendible to a 
larger time interval in this case too. □ 

There is another way of looking at this result. Suppose that (/, A,/i) 
is a regular solution of ( |1.1D -( [L5| ) on the whole space and that [0,T[ is its 
maximal interval of existence. Then by a similar argument to the above the 
solution extends in a manner to the set ([0,T] x IR^) \{(T,0)}. Thus if a 
solution of (|1.1|) -( |L5|) develops a singularity at all the first singularity must 
be at the center. 
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